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Abstract 

> ■ 

If^ • In this talk I give an overview of the work done during the last 15 

, years in collaboration with the late Adrian Patrascioiu. In this work 

^ I we accumulated evidence against the commonly accepted view that 

I ' theories with nonabelian symmetry - either two dimensional nonlinear 

. a models or four dimensional Yang-Mills theories - have the property 

I of asymptotic freedom (AF) usually ascribed to them. 

^ . 

> , 

: 1 Introduction 

> ■ 

^ I Our present view of nature is based on reductionism: bulk matter consists 

' of atoms, which consist of electrons and nuclei; nuclei consist of protons and 

neutrons, which in turn consist of quarks and gluons. This simple idea that 
something 'consists' of something smaller becomes increasingly inadequate 
as we go down; certainly the statement that the proton consists of three 
quarks has to be taken with a big grain of salt: depending on the energy 
applied, it also seems to consists of three quarks and an arbitrarily large 
number of quark-antiquark pairs as well as gluons. But most of all, all these 
'constituents' do not exist as particles in the usual sense: they cannot be 
isolated and do not fit Wigner's definition of a particle as the embodiment 
of an irreducible representation of the Poincare group. 
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The catchword to describe the last property is 'confinement' and the 
theory that is supposed to embody it is Quantum Chromodynamics (QCD), 
based on Yang-Mills theory, enriched with fermion (quark) fields. Lattice 
gauge theory at strong coupling shows indeed this property of confinement 
and a huge body of lattice simulations has given us confidence that QCD 
has a very good chance to describe correctly the spectrum of hadrons. Of 
course the lattice is an artefact and one would like to get rid of it by taking 
the so-called continuum limit. The way to do this is to search for a critical 
point of the lattice theory, in which the dynamically generated length scale 
diverges in lattice units. 

It is part of the conventional wisdom and stated in numerous textbooks 
that this critical point is located at vanishing (bare) coupling, and that the 
continuum limit hence enjoys a property called asymptotic freedom (AF), 
meaning that the effective interaction strength goes to zero with increasing 
energy. This prediction is based on computations done in perturbation theory 
(PT), but so far a proof does not exist. 

Similar predictions have been made about two-dimensional nonlinear a 
models with nonabelian symmetry group, which in many ways can be con- 
sidered as toy models for QCD. But it should be stressed that even for this 

toy version AF has neither been proved nor disproved mathematically - to 
do so remains an important challenge for mathematical physics. 

In fact Patrascioiu and the author have questioned the textbook wisdom 
and over the years we have accumulated a lot of evidence, both analytic 
and numerical, in support of a conjecture contradicting the conventional 
picture, namely that both four-dimensional Yang-Mills theory and its two- 
dimensional toy analogues, considered as lattice statistical systems, have a 
critical point at a finite value of the lattice coupling strength, which then 
corresponds to a non-Gaussian fixed point of the Renormalization Group 



In this lecture I will focus on the two-dimensional toy version, which 
is easier to study. Concretely I will consider 0{N) models (classical spin 
models) which are defined as follows: we consider configurations of classical 
spins s on a simple square lattice Z^, i.e. maps 



(RG). 
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for any finite subset A of the lattice 1? a Hamiltonian is given by 

= - s{^) ■ s{y) ; (2) 

On the configurations restricted to A, i^A induces a Gibbs measure 

^e-^^'-dMs}) . (3) 

As usual, one has to take the thermodynamic limit k Z' 1? . The mass gap 
is then defined as the inverse correlation length 

= = - lim ^ ln(s(0) ■ s[x)) . (4) 

The textbooks state that there is a fundamental difference between the 
models with N = 2 (abelian symmetry) and N > 2 (nonabelian symmetry), 
namely, whereas for iV = 2 we have 

m{(3) = for /3 > /3kt , (5) 

for > 2 on the contrary 

m{p) > for all > . (6) 

The first part of the statement, which goes back to the seminal paper of 
Kosterlitz and Thouless PP, has in fact been proven by Frohlich and Spencer 
[2j long ago. The second part, on the other hand, remains unproven to date 
and it represents an important open problem of mathematical physics to 
either prove and disprove it. Its importance for the understanding of two- 
dimensional ferromagnets is obvious, but maybe even more important is its 
analogy to the problem of mass generation in four dimensional Yang-Mills 
theory or QCD. QCD as a valid theory of strong interactions requires a mass 
gap due to the short range nature of the nuclear forces, and it needs a mass 
scale ('string tension') describing the strength of the confining force. 

2 Lattice construction of Quantum Field The- 
ories 

Let me briefly revisit the principles of constructing a massive (euclidean) 
Quantum Field Theory as a continuum limit of a lattice statistical model. 
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We assume that the thermodynamic hmit of a theory on the lattice Z has 
been taken aheady and yields an infinite volume translation invariant Gibbs 
state. In the high temperature (strong coupling) regime (3 << 1 there is 
exponential clustering as can be shown using convergent cluster expansions. 
So the system produces its own dynamical scale, the correlation length ^ 
defined in (jlj). 

We now change scale, using ^ as the standard of length; so the lattice 
appears now as the rescaled version of Z^, namely aZ^ with a = 1/C,- I 
want to stress that therefore the lattice spacing is not a freely chosen 
parameter of the model, but a dynamically determined quantity. 

In particular a continuum limit a ^ requires ^ oo, i.e. it requires the 
existence of a critical point. By taking the correlation length as the standard 
of length, one will produce a massive continuum limit (mass 1 by our choice). 

Returning now to the 2D 0{N) models, the procedure means that the 
continuum limit of the spin correlation functions is given by 

Sai,...aAxi, - ■ ■ ,Xn) = hm Z(/?) ""^^ (s^, (Xi^ , • • • , Sa„ (a^nO ) • 

where /3c denotes the critical value of (3 (finite or infinite) and Z{(3) is a field 
strength renormalization that has to go to for (3 —>■ (3^-, if the continuum limit 
is to be a nontrivial Quantum Field Theory - because in the continuum limit 
the correlation functions (Schwinger functions) has to blow up at coinciding 
points. A simple possible choice for Z[(5) is 

m = I , (8) 

(cf. for instance j3j) where x is the susceptibility of the spin model: 

x = Y.m-s{^))- (9) 

While the construction of the continuum limit works whether (3c is finite 
or infinite, a finite value for the 0(3) model has the crucial consequence that 
it is not aymptotically free: 

/? - c < oo ^ No AF ! (10) 

This follows from the bounds in 4] and as has been stressed in [H]. 
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It should be noted that the construction of the continuum hmit in AD 
Yang-Mills theory follows the same principle; instead of the correlation length 
^ one may take 1 / ^/a where a is the string tension giving the slope of the 
confining potential - it is believed that these two scales are equivalent, i.e. 
their ratio tends to a finite value as the continuum is approached. In any 
case, with this choice of scale (assuming it exists) one obtains a confining 
continuum limit. In this sense one may say that the 'confinement problem' 
has been solved by lattice gauge theory. 

In all these cases {2D nonlinear a models and AD gauge theories) the con- 
tinuum limit has no free parameter corresponding to a choice of the coupling 
constant, there is only a scale (multiple of C,) that may be chosen. This is the 
famous 'dimensional transmutation' that is normally ascribed to properties 
of the perturbation theory for these models; we see that it is a much simpler 
and more general fact. 



3 Arguments for Asymptotic Freedom 

All arguments are based on perturbation theory (PT). This is a formal low 
temperature expansion in powers of 1/(3 = g. For the 0{N) models it predicts 
the flow of the running coupling to be governed by the Callan-Symanzik (3 
function 

feW = -^-^+0(«'), (11) 

see [0]). 

PT is in principle just the application of Laplace's method for the deriva- 
tion of asymptotic expansions for sharply peaked integrands. But one should 
be wary because we are dealing with infinite systems and the usual theorems 
(see for instance [3 El) do not apply. The idea is the following: the Gibbs 
factor exp{—(3H/^) for /? >> 1 will be sharply peaked around the ground 
state configuration of H\. The ground state of an 0{N) model is simply any 
configuration in which all spins ar equal: 

tttt ••• TT (12) 

The expansion is produced by treating the configurations as small fluctua- 
tions around such a ground state. For a finite system (L < oo) this procedure 
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is clearly correct; but in an infinite system (L = oo) there will always be large 
fluctuations according to the Mermin- Wagner theorem jHl QHl • So this 
means that PT should be considered as a priori dubious. 

Traditionally one has simply done the expansion in a finite volume (or 
with some other infrared regulator) and then considered the limit of the 
expansion coefficients as termwise as the the thermodynamic limit is taken (or 
the regulator is removed). In 1980 Bricmont et al ^2] actually showed that 
for invariant observables in the 0(2) model, this formal procedure produces 
a valid asymptotic expansion for expectation values in the infinite volume. 
Nothing of that sort has been achieved for N > 2. 

But confidence was boosted by the finding of Elitzur and later David 
that for N > 2 the termwise limit exists. This is generally (but with- 
out justification) assumed to mean that one obtains indeed the correct low 
temperature asymptotic expansion by PT. 

But the problem cannot be settled by looking at the termwise limits. 
The structure of the problem is as follows: let A be some observable; then 
for L < oo we have 

M 

{A)l = 5^c„(L)r" + ^M(/3,^) (13) 

n=0 

with 

\Rm{(3,L)\<Km{L)[3-^'-\ (14) 

Elitzur and David showed that lim^^^oo c„(L) exists, but to justify PT one 
would have to prove uniform bounds on Km{L). 

4 Reasons for doubt 

I list a few facts that should be reason not to trust PT in the 0{N) models 
for > 3, 

• I mentioned already that the Mermin- Wagner theorem requires large 
fiuctuations at all temperatures and so destroys the logical basis of PT. 

The intuitive reason for this fact lies in the existence of localized defects 
of arbitrarily low energy, which disorder the system at all temperatures. 
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We dubbed these defects 'superinstantons', because they will 'beat' the 
instantons (which have fixed energy) at low temperature. They can 
be simply described as follows: fix a spin at the origin to a certain 
direction, say e/v, the unit vector in the A^th direction; also fix all spins 
at approximate distance R from the origin to another direction e*. The 
configuration of minimal energy with these boundary conditions (be) 
is called superinstanton of size R (and rotation angle arccos(e/v ■ e*)). 
It turns out that the energy of such a superinstanton is 

E{R) = 0{R~^) in D = 1 

E{R) = 0(1/ \nR) mD = 2 

E{R) = 0(i?°) mD>2 (15) 

(see for instance [IHl)- A notable property of these excitations is that 
they do not exploit the fact that the spin has components; they are 
0(2) like for all N and all D. They are local minima of the energy 
(with the prescribed be); in the continuum the 0(2) superinstanton 
can be seen to be the conjugate harmonic function to the vortex that 
plays such a fundamental role in the Kosterlitz-Thouless (KT) theory. 

• Ambiguity of PT [IH]: 

The thermodynamic limit of 2D 0{N) models ought to be independent 
of be. This is plausible because of the Mermin- Wagner theorem, but it 
is a definite fact in the massive phase (which for N > 2 according to the 
conventional wisdom is all there is) because of exponential clustering. 

Nevertheless Patrascioiu and the author found some be that lead to PT 
coefficients that differ from the standard ones (obtained with periodic 
be) in the termwise thermodynamic limit. These be were inspired by 
the superinstantons and called superinstanton be (sibc). They are de- 
fined essentially by enforcing a superinstanton of rotation angle zero. If 
we consider 0{N) invariant observables, sibc can be defined by fixing 
the spin at the origin and all the spins at the boundary of a square to 
a certain fixed direction. 

It is inportant to realize that these are indeed legitimate be for the 
model: in the thermodynamic limit the effect of the boundary spins 
disappears, and fixing one spin in the center is inconsequential, since 
we are considering invariant observables. 
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For N > 2 (but not for = 2!) PT at one-loop level (termise thermo- 
dynamic limit), already produces an answer different from the standard 
one if one uses sibc . This means that at least one of the two results, 
possibly both, are incorrect, since the true (nonperturbative) thermo- 
dynamic limit is the same and an asymptotic expansion is unique, if it 
exists. 

• Richard's truncated sphere model ^7]: 

If one modifies the 0(3) model by restricting the z component of the 
spin to satisfy the constraint < 1 — e, the system becomes massless 
at sufficiently low temperatures. This was shown by J.L.Richard by 
the use of Ginibre's correlation inequalities for the x, y components 
and comparing the model to an 0(2) model. 

But PT would tell a different (and untrue) story: PT, expanding around 
a ground state oriented towards the equator does not see the constraint 
at all. It yields the Callan-Symanzik (3 function of the ordinary 0(3) 
model (fTTjl and therefore, by the usual reasoning, says that the model 
is aymptotically free and massive for all (3\ 

5 The case against AF: discrete and continu- 
ous models 

It is interesting to make a 'lattice approximation of the target space' by 
replacing for instance the 2-sphere 5*2 of the 0(3) model with (the set of 
vertices of) a platonic solid such as the icosahedron or dodecahedron. In 
both of these cases the symmetry is also reduced from the continous 0(3) to 
the discrete icosahedral subgoup Y. 

The study of the dodecahedron model was begun already long ago, both 
analytically and numerically in UHJ and continued in a detailed 
numerical study of the icosahedron model was carried out in [3]. The picture 
that emerged is that the dodecahedron has a critical interval between its 
massive high and low temperature phases (similar to the abelian Z„ models 
for large enough n), whereas the icosahedron model has an isolated critical 
point Pc- In both cases our main interest is in the (massive) continuum limit 
arising by approaching criticality from the high temperature phase. 
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Numerically we found evidence for symmetry enhancement from Y to 
0(3) in both cases, suggesting that the massive continuum limits of the 
discrete models lie in the same universality class as the 0(3) model. 

Concretely, we considered the following objects: 



Renormalized spin-spin correlation in momentum space 



where the data presented in ^201 suggest 



(16) 



lim G,.(p)|dodeca = lim Gr{p)\o{3) 

5— >oo g-^oo 



(17) 



Renormalized spin-spin correlation in x space 



where the data presented in |2I] suggest 



(18) 



lim Gr{x) 

5— >oo 



ico = lim G^(a;)|o(3) 

5— ►CXD 



(19) 



Renormalized 4-point coupling constant 



where the data in [21] give evidence for 



(20) 



(21) 



This suggests the overall conclusion that the dodecahedron, the icosahe- 
dron and the 0(3) model all have the same massive continuum limit. This is 
relevant for our main issue, because the discrete models cannot be expected 
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to be asymptotically free. In fact the running coupling proposed by Liisclier, 
Weisz and Wolff [221 



in the icosahedron model has found to have an ultraviolet fixed point 

595 ^0. (23) 

So assuming the universality suggested by our does indeed hold, this 
means that the 0(3) model cannot have AF. 



6 The case against asymptotic freedom: per- 
colation 

The strongest argument against AF developed by Patrascioiu and the author 
is based on the analysis of certain percolation properties. The idea was 
presented first already in 1991 and developed in a number of papers [2211211 
I2nil2ni- Since the argument does not constitute a rigorous proof, a detailed 
numerical study to bolster it was carried out much later [S]. 

The argument might seem tricky, but it involves some solid analytic re- 
sults, on which the heuristic/numerical arguments builds. 

• The first step are two modifications of the model that should not change 
its universality class: 

— We replace the square lattice 1? by a triangular lattice T; this is 
achieved simply by adding extra bonds along one of the diagonals 
of each elementary square. All the bonds carry the same Gibbs 
factor, so the model is really living on an isotropic triangular lat- 
tice. There is no question that this modification does not affect 
the continuum limit, i.e. does not change the universality class. 

— Next we introduce a constraint in the Gibbs measure that limits 
the angle between neighboring spins to some maximal value: 

e-^^ ^ e-^^ JJ ^(s(x) • s(y) - c) 

{xy) 

('standard action') — > ('cut action') (24) 
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For large values of (3 this change is completely innocuous, because 
the standard action Gibbs factor will already make large angles 
extremely unlikely. Our arguments will be completely indepen- 
dent of the value of /3, so one might assume that it is already fixed 
at some large value. It turns out, however, that one can change 
the point of view by putting /3 = and varying the cut parameter 
c instead, since the presence of the constraint also has a ferro- 
magnetic ordering effect of the system. One finds indeed that the 
correlation length increases with increasing c and a critical value 
of c somewhere near 0.76 seems to exist. 

We also performed a direct numerical test of universality between 
this 'cut' model with /5 = and the standard action model by mea- 
suring the so-called 'step scaling function' that gives the change 
of the LWW coupling under doubling the scale as a function of 
the LWW coupling itself. The data show quite good agreement 
between the two models (see 0). 

• An Ising model is imbedded in the 0(3) model by setting 

(T^ = sg\is^{x) (25) 

as is done for the well-known cluster algorithms. Note that the def- 
inition of the Ising spins is well-defined except on a set of measure 
zero. 

• A correlated bond percolation model is set up following the original 
work of Fortuin and Kasteleyn |27J and adapted to 0{N) models by 
Wolff p8j: 

Given a spin configuration, bonds (nearest neighbor pairs) are activated 
with a conditional probability p{{xy)\{s{x) x G T} that is determined 
by the change in the Gibbs factor when the z component of one of the 
spins is reflected: 

p{{xy)\{s{x) X G T} = e{a,ay) [l - e(/5«»-(«-^"'(3/)-«-(3/)))] (26) 

where 

s RzS = s — 2ez{ez ■ s) (27) 

After activating the (xy) bonds independently with this probablity one 
forms connected bond clusters; then one averages over all spin config- 
urations with the Gibbs measure. 
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A theorem that is essentially due to Fortuin and Kasteleyn j27j now 
relates the expected size {Cfk) (i-e. number of vertices) of the bond 
cluster attached to the origin to the susceptibility of the imbedded Ising 
spins: 

Xising = {Cfk) ■ (28) 
This theorem has an important corollary: 



If (Cfk) = c)0,then C, = oo 



(29) 



• The next step is to divide the 2-sphere into 3 regions: 

— Equatorial strip = {s G 5*^1 \sz\ < |} 

— North polar cap P+ = {s G S'^ls^ > f} 

— South polar cap = {s G S'^\sz < — |} 

The idea is now to study clusters of S^, U V~ = V^- 

The activation probability for the bond {xy) is always 1 if c > 1 — y 
and s{x),s{y) G (or s{x),s{y) G V~), because in this case flipping 
one of the spins would violate the constraint. So we have 

(Cfk) > (Vt) (30) 



• The main result of [221 12^] is the following: 

If for some e and some c > 1 — e^/2 does NOT percolate, 
then ^ = oo. 

The idea of the proof in intuitive terms is as follows: assume that cS'^ 
does not percolate, then there are two possibilities: 

— percolates =^ V^, V~ percolate: but this is impossible because 
in 2D there cannot be two disjoint percolating clusters (actually 
this is a 'principle', proven only for Bernoulli percolation and the 
Ising model) 

— Vt does not percolate, but prevents percolation of 

=^ V^, form 'rings' of arbitrarily large size; neither nor 
its complement percolate. Now by a lemma of Russo j^Hl, if on 
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a self-matching lattice (such as T) neither a set nor its comple- 
ment percolate, the expected cluster size of both of them diverges: 
{V'^) = 00^ {Crx) = 00 ^ ^ = 00, which completes the argu- 
ment. 

• So the question remains whether it is possible for an equatorial strip 

2 

to percolate for arbitrarily small e and c > 1 — y- This is a priori hard 
to imagine. In any case, in [221 12nj we gave an argument that leads to 
masslessness even if we assume this implausible situation to occur: 

Assume percolates. Then the same is of course true for any S^' with 
e' > e. Taking e' close to 2 one obtains an 'ocean' with only rare and 
small islands corresponding to V'^. Focussing on the two components 
Sx, Sy, and taking (3 very large, we have a low temperature 0(2) model 
with fluctuating coupling on that ocean. One can arrange for the is- 
lands to cover an arbitrarily small fraction of the lattice and at the 
same time make the effective temperature for the 0(2) model arbitrar- 
ily small. So one expects that one ends up in a massless Kosterlitz- 
Thouless (KT) phase. The only problem is that the Frohlich- Spencer 
proof has not been adapted for this situation. 

Of course I should stress that I don't think this situation ever arises; 
I rather expect that at low enough temperature (large enough c) the 
clusters of an arbitrarily small polar cap form rings of any size and have 
divergent cluster size. This conforms to the superinstanton picture and 
parallels the situation rigorously established for the 'cut' 0(2) model 
(see jMEOl)- 

• Finally I want to briefly describe the extensive numerical study of the 
percolation properties in 0: 

For simplicity we studied only the model with (3 = 0, but varying cut 
c. A nonzero (3 would only order the system more; if it is massless at 
P = it would a fortiori be so for /5 > 0. We scanned the (c, e) plane 
and measured the ratio r(L) = {V^)/{S^) as a function of the lattice 
size L. For fixed 0.78 < e < 1 and varying c betwqeen -0.1 and 0.9 
three regions could be distinguished: 

— For small c, r(L) increases rapidly with L, indicating percolation 
of - it should be noted that in this regime the clusters of V'^ 
and V~ are allowed to touch. 
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— For intermediate c, r{L) decreases rapidly with L, presumably 
converging to and indicating percolation of (Sg. The boundary 
between these two regions is a rather sharply defined value of c, 
where r(L) is practically independent of L 

— for even larger c, this L dependence of r(L) seems to disappear, 
or even change direction again, but at most showing a very mild 
increase with L (maybe powerlike). We take this as an indication 
that now neither of the sets Ve, percolates (which means by 
Russo's lemma that they both have divergent mean cluster size). 

It turns out furthermore that for c < 0.76 or so the intermediate region 
of percolating seems to disappear altogether. To assess the credibil- 
ity of these conclusions, one should note that we went to quite large 
lattices (up to L = 1280), but numerical results are of course never a 
substitute for a proof. 

We can combine the picture suggested by the numerical results into a 
(semiquantitative) 'percolation phase diagram' (Fig. taken from [S] . 

To corroborate our conclusion that for c < 0.76 neither Ve nor 
percolates, we carried out a different analysis of the equatorial clusters 
with width e = 0.75; this time we compared them to the clusters of 
V^, chosen such that the polar cap has the same area as the equatorial 
strip iS^; also both sets cover the same fraction of the lattice. The ratio 
r'(L) = (V^) / {S^) is always larger than 1, presumably because the 
polar cap has less boundary than the equatorial strip of the same area 
(in agreement with a conjecture stated in j3I])- For small c (roughly 
up to c = 0.4) r'(L) shows a rapid increase with L, indicating that Vd 
forms rings (of all sizes) whereas the clusters of have finite mean size. 
At c ~ 0.4 the behavior of r'{L) changes: it is still growing with L, but 
now only powerlike. This can only indicate that both kinds of clusters 
now form rings of arbitrary size; it rules out for all practical purposes 
that Se percolates, which would require a decrease to as L ^ cxd. 

The numerical results of this study are shown in Fig|^ taken from [S]. 

To sum up: our numerical study gives strong evidence that for suffi- 
cently small e (about < 0.76) does not percolate for any value of c. 
This implies according to our results that for c > 1 — 0.76^/2 = 0.745 
the correlation length ^ = oo. The real challenge is of course to give a 
mathematical proof of this. 
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-1 -0.5 0.5 1 

c 

Figure 1: Percolation phase diagram of the 0(3) model on the T lattice. C 
is the line c = 1 — above the dashed line D percolates, above the 
dotted line (S^) < oo. 

7 Where do we stand? 

Clearly, the problem of AF is still open 

• Textbook wisdom is insufficient to settle it positively 

• Counterarguments are not rigorous 

• Any progress towards proving or disproving AF is very desirable. For 
this reason programs such the one of K.R.Ito attempting to prove 
AF by a controlled Renormalization Group approach, are very welcome. 

A final remark: Solving the problem of AF for the 2D toy models which 
were mostly discussed here would be a big step towards understanding it 
also in 4D Yang-Mills theory and thereby also towards a solution of one of 
the Million Dollar Problems posed by the Clay Mathematics Institute. 
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Figure 2: The ratio of the mean cluster size of a polar cap of height 0.75 to 
that of an equatorial strip of the same height 
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